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We present theoretical and dynamic Monte Carlo simulation results for the mobility and micro- 
scopic structure of 1+1-dimensional Ising interfaces moving far from equilibrium in an applied field 
under a single-spin-flip "soft" stochastic dynamic. The soft dynamic is characterized by the prop- 
erty that the effects of changes in field energy and interaction energy factorize in the transition rate, 
in contrast to the nonfactorizing nature of the traditional Glauber and Metropolis rates ("hard" 
dynamics). This work extends our previous studies of the Ising model with a hard dynamic and 
the unrestricted SOS model with soft and hard dynamics. [P. A. Rikvold and M. Kolesik, J. Stat. 
Phys. 100, 377 (2000); J. Phys. A 35, L117 (2002); Phys. Rev. E 66, 066116 (2002).] The Ising 
model with soft dynamics is found to have closely similar properties to the SOS model with the 
same dynamic. In particular, the local interface width does not diverge with increasing field, as it 
does for hard dynamics. The skewness of the interface at nonzero field is very weak and has the 
opposite sign of that obtained with hard dynamics. 
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I. INTRODUCTION 



The structure and dynamics of surfaces and interfaces significantly influence a host of material properties. Conse- 
quently, an enormous amount of work has been devoted to the study of moving and growing interfaces 1, ^. However, 
despite the fact that many important interface properties, such as mobility and chemical activity, are largely deter- 
mined by the microscopic interface structure, the bulk of this effort has concentrated on large-scale scaling properties. 

Since the detailed physical mechanisms of the interface motion are most often unknown, it is useful to model the 
dynamic as a stochastic process defined by a set of transition probabilities. It is therefore important to gain better 
insight into how the driving force (such as an applied magnetic or electric field or a chemical-potential difference) may 
alter the microscopic interface structure for different stochastic dynamics. Recently we have studied the influence 
of the stochastic dynamics on the microscopic structure and mobility of Ising and solid-on-solid (SOS) interfaces 

ODQ 



that move under two types of Glauber dynamics 
Hamiltonian 



Both Ising and SOS interfaces are described by the Ising 



H = — S x>y (J x S x +l,y + JyS x ,y+l + H) , (1) 

where s XjV = ±1 is an Ising spin at lattice site (x,y), ^2 X runs over all sites on a square lattice, and J x and J y 
are ferromagnetic interactions in the x- and y-directions, respectively. The quantity H is the applied "field," and the 
interface is introduced by fixing s x . y = +1 and —1 for large negative and positive y, respectively. We take H > 0, 
such that the interface on average moves in the positive y direction. (This model is equivalent to a lattice-gas model 
with local occupation variables c x>y € {0, 1}, see details in Ref. |5j.) The difference between the two interface types is 
that the Ising interface allows overhangs and bubbles, while these are forbidden in the SOS interface. However, at low 
and intermediate temperatures overhangs and bubbles in an Ising interface are rare, and a short interface segment is 
likely to be indistinguishable from an SOS interface. A typical SOS interface is illustrated in Fig. 

In addition to bubbles that are generated at the interface by pinching off of protrusions or indentations, Ising 
models in an applied field can also contain bubbles created by homogeneous 

While such bubbles destroy the integrity of the interface at very strong fields, they have only a minor influence on 
the mobility at moderate fields In this paper, like in Ref. (3, we exclude such nucleated bubbles by setting the 

transition rate equal to zero for sites that have no neighbors with the opposite spin direction. As a consequence, the 
bulk phases far from the interface are uniform. 

The interface dynamic is defined by the set of single-spin transition probabilities, W(s Xiy — > —s Xi y), and time is 



2 



measured in units of attempted MC updates per site (MC steps per site, or MCSS). The first one of the dynamics used 

^ is the standard discrete-time Glauber dynamic with the transition probability [j]] 



in the aforementioned studies 



%(««,»-»-«.,»)= [l + e^ E ] 1 , (2) 



where AE is the total energy change that would result from the transition. Although AE can be written as a sum of 
the energy change AEj, due to the change in the interaction part of the Hamiltonian, and AEh, due to the change in 
the field energy, this transition rate itself cannot be factorized into a product of parts that depend only on AE.j and 
AEh, respectively. This dynamic is therefore classified as "hard" in the literature on driven particle systems 

The second type of dynamics is defined by transition probabilities that factorize into a part that depends only on 



Ej a 



AEj and one that depends only on AEh- Such dynamics are known as "soft" 10] . In our recent study of a driven 
SOS interface with soft dynamics we used, for reasons of mathematical convenience, the "soft Glauber dynamic" 
in which each of the two parts has the Glauber form: 

W SG (s XtV - -a WlV ) = [1 + e" AE -} - 1 ■ [1 + e^] 1 . (3) 

Soft dynamics (usually with the field part proportional to a Metropolis transition rate ^) are often used for 
lattice-gas simulations, in which the field term corresponds to the entropic part of a chemical-potential difference. In 



Ref. 4| we showed, in agreement with a theoretical prediction in Ref. |3j, that the soft dynamic leads to a microscopic 
SOS interface structure that is identical to the equilibrium interface in zero field, irrespective of the value of the 

n 

applied field. This is in contrast with the result for hard dynamics |5J, which lead to an intrinsic interface width that 
increases dramatically with the field. The purpose of the present paper is to study the effects of the soft Glauber 



dynamic for an Ising interface (which may contain bubbles and overhangs) and compare those with the cases of an 
Ising interface with the standard (hard) Glauber dynamic [j| and SOS interfaces with the soft Glauber dynamic 

n 

and hard Glauber dynamic |5|. 

The rest of this paper is organized as follows. Theoretical results for the interface structure and velocity are surveyed 
in Sec. ^ Comparisons with extensive dynamic MC simulations are given in Sec. 1 1 1 1 1 with results for the interface 
velocity in Sec. IIII Al and for the interface structure and skewness in Sec. lIIIBl Our conclusions are drawn in Sec. lIVI 



II. INTERFACE STRUCTURE AND VELOCITY 



With the Ising Hamiltonian there is only a finite number of different values of AE, and the spins can therefore be 
divided into classes [13! Ljl ll| , labeled by the spin value s and the number of broken bonds between the spin and its 
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nearest neighbors in the x and y direction, j and k, respectively. The eighteen different fsing spin classes are denoted 



jks with j, k £ {0, 1 , 2}. They are listed in Tabled and subsets are also listed in Table ITT1 and shown in Fig. ^ 



The Burton-Cabrera- Frank SOS model [lfj considers an interface in a lattice gas or 5=1/2 Fsing system on 
a square lattice of unit lattice constant as a single-valued integer function h(x) of the cc-coordinate, with steps 
S(x) = h(x + 1/2) — h(x — 1/2) at integer values of x. A typical SOS interface configuration is shown in Fig. ^ 
The heights of the individual steps are assumed to be statistically independent and, in the case of a flat interface, 
identically distributed. These assumptions are exact in equilibrium |16| . The step-height probability density function 
(PDF) is given by the interaction energy corresponding to the \S(x)\ broken J^-bonds between spins in the columns 
centered at (x — 1/2) and (x + 1/2) as 



p[6(x)] = Z^y 1 ^^ e^ 6 ^ 



(4) 



The Boltzmann factor X = e~ 2/3J:r determines the width of the PDF, and the Lagrange multiplier 7(0) maintains the 
mean step height at an ^-independent value, (S(x)) = tan <fr, where <fi is the overall angle between the interface and 
the x axis. The Lagrange multiplier is given by 



.7(0) _ 



(1 + X 2 )tan0+ (l - X 2 ) 2 tan 2 $ + AX'< 



1/2 



2X (1 + tan, 



The partition function for the step height S(x) is 



Z(cb) = £ X 



1-X< 



(See details in Refs. 



11) Si 



5= — oo 



1 - 2X cosh 7 (0) +X 2 



mple results are obtained for = 0, which yields 7(0) = and 



(5) 



(6) 



Z{Q) = {1+X)/{1-X) 



(7) 



and for = 45°, which yields 



^ 45 °)=(1+1 2 )/2X 



(8) 



and 



Z(45°) = 2(l+X 2 )/(l- X 2 ) 



(9) 



For soft dynamics (but not for hard dynamics), X remains independent of H when the system is driven away from 
ill 



equilibrium 
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The mean spin-class populations, (n(jks)}, are all obtained from the product of the independent PDFs for 5{x) and 
S(x+1). Symmetry of under the transformation [x, <ft, S) — > (— x, — <f>, — 6) ensures that (n(j'fe— )) = (n{jk+)) 

for all j and k. Numerical results illustrating the breakdown of this up/down symmetry for large \H\ are discussed 
in Sec. IIII Bl As discussed in Ref. t 3j, calculation of the individual class populations is straightforward but somewhat 
tedious, especially for nonzero <j>. The final results are summarized in Table UTI 

Whenever a spin flips from —1 to +1, the corresponding column of the interface advances by one lattice constant 
in the y direction. Conversely, the column recedes by one lattice constant when a spin flips from +1 to — 1. The 
corresponding energy changes are given in the third and fourth columns of Table [I] Since the spin-class populations 
on both sides of the interface are equal in this approximation, the contribution to the mean velocity in the y direction 
from sites in the classes jk— and jk+ becomes 



(v y (jk)) = W {0AE(jk-)) - W (/3AE{jk+)) 



(10) 



The results corresponding to the soft Glauber transition probabilities used here, Eq. (|3J), are listed in the last column 
of Table ITU The mean propagation velocity perpendicular to the interface becomes 



(v ± (T, H, </>)) = cos cf> X>C/fcs)) H W) 

3,k 



(11) 



where the sum runs over the classes included in Table [Tj] While the general result is cumbersome if written out in 
detail, using the fact that e~ 4/3Jl = X 2 for the soft Glauber dynamic [], J, J, we obtain relatively compact formulas 
for the special cases of <f> — and = 45° : 



(v ± (T,H,0)) = X 



1 



X 



l + X 2 {1 + X) 2 (1-X 2 ) 



X 2 



2(1 + 2X) 

1 + e^Jy + 1 + X 2 e i ' 3J y 



tanh(/?ff) 



(12) 



and 



(vi(T,ff,45°)) = 



2X 2 



(l + X 2 ) 2 1-X 4 



1 + 2X 2 + 3X 4 



1 + e 



4/3 Jy 



X 4 



1 + X 2 e 4 P J y 



tanhC&ff) 
2V2 



(13) 



III. COMPARISON WITH MONTE CARLO SIMULATIONS 



We have compared the analytical estimates of the propagation velocities and spin-class populations developed above 
with MC simulations of the same model for J x = J y = J. The details of our implementation of the discrete-time 

5 



Q are described in Ref. ^] 



By keeping only the interface sites in memory, 



n-fold way rejection-free MC algorithm |13j are described in Ref. 
the algorithm is not subject to size limitations in the y direction, enabling simulations for arbitrarily long times. 

Our numerical results are based on simulations with L x = 10 000 and <p between and 45° for several temperatures 
below T c . (T c — — 2J/ln(v / 2 — 1) w 2. 269 J is the critical temperature for the isotropic, square-lattice Ising model 
[l7|.) In order to ensure stationarity we ran the simulation for 10 000 n-fold way updates per updatable spin (UPS) 
before taking any measurements, and the results were averaged over 200 000 UPS 18j . 

A. Interface velocities 

First we compare the simulated interface velocities with the analytical approximation, Eq. (flip . Figure shows the 
normal velocity vs H for <j) — and a range of temperatures up to T c . There is excellent agreement between the MC 
results and the theory for temperatures below 0.8T C . 

The dependence of the normal velocity on the tilt angle cf) is shown in Fig. [3] for several values of H/J between 
0.1 and 3.0. At T — 0.2T C the velocity increases with cj> m agreement with the polynuclear growth model P, 2, 2^| 
at small angles and the single-step model for larger angles 

D Q Q [Fig. Eta)]. At T = 0.6T C , on the other 
hand, the velocity is nearly isotropic, with a weak increase with <p for the strongest fields [Fig.|3Jb)]. For the lowest 
temperature the agreement between the simulations and the analytical results is excellent everywhere. For the higher 
temperature it is also reasonable, but better for weak than for strong fields. 

The temperature dependence of the normal interface velocity is shown in Fig. 0] for several values of H/J between 
0.1 and 3.0. The agreement between the simulations and the analytical results is excellent except for combinations 
of high temperatures and strong fields. In contrast to the results for hard dynamics (see Fig. 5 of Ref. P for Ising 

n 

interfaces and Fig. 8 of Ref. 5] for SOS interfaces), the velocity goes to zero at T = for all values of H, not just 



for H/J < 2. This result agrees with our finding for the SOS model with soft dynamics |4|. As predicted by the 
theoretical results in Ref. Q , there is thus no discontinuity in the interface velocity at T — and H/ J = 2 for soft 
dynamics. 

B. Spin-class populations and skewness 

A closer look at the performance of the mean-field approximation for the interface structure is provided by the 
mean spin-class populations. The analytical predictions for the class populations are based on the assumption that 
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different steps are statistically independent. A comparison of the simulation results with the analytical predictions 
therefore gives a way of testing this assumption. 

The ten mean class populations that have nonzero populations in the SOS approximation, (n(01s)), (n(lls)), 
(n(10s)), and (n(21s)), and (ra(20s)) with s = ±1 are shown vs H in Fig. EJa) for tj> = and T = 0.6T C . Filled 
symbols represent s = +1, while empty symbols (almost completely hidden behind the corresponding filled symbols) 
represent s = — 1. The class populations are practically independent of H, in agreement with the theoretical prediction 



for the soft dynamic 



i 



and in contrast to the result for the Ising model with hard Glauber dynamics (see Fig. 7(a) 
of Ref. |3j). Deviations from the SOS approximation are indicated by the nonzero populations in the classes with two 
broken y-bonds, (n(12s)), (n(22s)), and (n(02s)), which are shown in Fig. [SJb). These populations are only of the 
order of 10~ 3 , about two orders of magnitude less than for the hard Glauber dynamic (see Fig. 7(b) of Ref. [3j), and 
they show significant differences between the two spin values (see below). Figure|SJc) shows the combined populations 
in classes with one, two, three, and four broken bonds, respectively. The results are dominated by the SOS-compatible 
classes and show good agreement between simulations and theory. 

The skewness between the spin populations on the leading and trailing edges of the interface are a consequence 
of short-range correlations between neighboring steps, and it is quite commonly observed in driven interfaces. This 
is the case, even when the long-range correlations vanish as they do for interfaces in the KPZ dynamic universality 
, to which the present model belongs. Skewness has also been observed in several other SOS-type models 
The correlations associated with the skewness generally lead to a broadening of protrusions on the leading 



class L 


L|23 


24 


25 





edg e ("hilltops"), while those on the trailing edge ("valley bottoms") are sharpened [24j, or the other way around 
[26j . In terms of spin-class populations, the former corresponds to (n(21— )) > (n(21+)) and (n(ll+)) > (n(ll — )). 
The relative skewness can therefore be quantified by the two functions, 

_ (n(21-)) - (n(21+)) 



introduced by Neergaard and den Nijs 



' (n(21-)) + (n(21+)) ' 
and 

_ (n(ll+)> - <n(ll-)) 



(14) 



(15) 



(n(ll+)> + <n(ll-)> ' 

These two skewness parameters for the current system are shown together versus H in Fig.[|Ja). The skewness is very 
weak, but it is not zero in contrast to the SOS model with soft Glauber dynamic Both these skewness parameters 
have the opposite sign and are about two orders of magnitude smaller than in the Ising model with the hard Glauber 
dynamic 3] [Fig.[(Jb)] and in the SOS model with the same hard dynamic [5j [Fig. Efc)]. 
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The skewness parameters p and e depend on spin classes which have nonzero populations in the SOS picture, and 
they can therefore be applied to both Ising and SOS interfaces. In the Ising case, however, much more pronounced 
differences are seen in in the populations of those classes that are not populated in the SOS model [see Fig.[5£b)]. For 
Ising models, these classes can be used to define further skewness parameters, such as 

(n(22-)) - (n(22+)) 

(n(22-)) + (n(22+)) ' 1 ' 

Here class 22— represents isolated bubbles of the metastable phase that persist as a "wake" behind the moving 
interface, while 22+ corresponds to a "bow wave" of bubbles of the stable phase in front of the interface, which are 
created by pinching-off of protrusions. Although the total density of such bubbles is about two orders of magnitude 

Q 

smaller with the soft dynamic studied here than with hard dynamics (compare Fig. Efb) with Fig. 7(b) of Ref. |3|]), 
the relative asymmetry parameters k are comparable, as shown in Fig. [J] Not surprisingly, for strong fields the entire 
bubble population is found in the wake, yielding k < 1 in this limit. 



IV. CONCLUSION 



In this paper we have continued our study of the dependence of the local structure of driven interfaces on the applied 



say 



The local interface 



field and temperature and on the form of the stochastic dynamics under which they move 
structure is of interest because it is it, rather than the large-scale scaling behavior, which determines such important 
interface properties as mobility and chemical reactivity. In particular, we studied the differences between interfaces 
moving under soft stochastic dynamics, in which the influences of changes in the field energy and the interaction energy 
factorize in the transition probabilities, and hard dynamics, which do not possess such a factorization property. 
We find that the results for the Ising model with the soft Glauber dynamic, which is the main topic of this study, 



differ relatively little from the SOS model with soft dynamics studied in Ref. 



Q. 

ing 1 and SOS 



In particular, the local interface 



width does not diverge with increasing H , as it does for both the Ising jy] and SOS 5] models with hard dynamics. 
As a result, the soft dynamics do not produce the discontinuity in the interface velocity at H/ J = 2 and T = that is 
seen for hard dynamics. The main qualitative difference between the Ising and SOS models with soft dynamics is that 
the interface skewness for the Ising model in nonzero field is not exactly zero, as it is for the SOS model. However, 
the skewness parameters p and e, which are based on SOS-compatible spin classes, are about two orders of magnitude 
smaller and have the opposite sign than what is seen for the hard dynamics. In contrast, the relative asymmetry in 
the populations of bubbles behind and in front of the interface (which do not occur in SOS models) can be expressed 
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ICO 

5], 



by the skewness parameter k and is comparable for the two dynamics. However, the absolute bubble density is about 
two orders of magnitude smaller with the soft than with the hard dynamics. Although a successful mean-field theory 
for the interface mobility of different models and under different stochastic dynamics was developed in Refs. 
comparable theory that predicts the skewness is still not available. 

Two important conclusions can be drawn from our studies. First, the strong differences between hard and soft 
dynamics make it evident that great care must be used in formulating and interpreting stochastic models of dynamic 
systems. Second, experimental observation of the field and temperature dependences of interface mobility and local 
interface structure could contribute significantly to devising correct stochastic models of nonequilibrium physical 
phenomena. 
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TABLE I: The spin classes in the anisotropic square-lattice Ising model. First column: the class labels, jks. Second column: 
the total energy per spin, E(jks), relative to the state with all spins parallel and H — 0, Eo = —2(J X + J y ). Third column: 
the change in the field energy resulting from spin reversal from s to — s, AEh {jks). Fourth column: the corresponding change 
in the interaction energy, AEj(jks). In columns two and three, the upper (lower) sign corresponds to s = —1 (s = +1). The 
first three classes have nonzero populations in the SOS approximation, and flipping a spin in any of them preserves the SOS 
configuration. The two classes marked f also have nonzero populations in the SOS approximation, but flipping a spin in any 
of them may produce an overhang or a bubble. The two classes marked X are not populated in the SOS approximation, but 
flipping a spin in any of them may produce an SOS configuration. The two classes marked * correspond to a bulk spin that 
is either parallel or antiparallel to all its neighbors. Flipping a spin in class 22s yields a spin in class 00— s. The transition 
probabilities for all classes except 00s (from which transitions are forbidden with the dynamic used here) are given by Eq. 10 . 



Class, jks 


E(jks) - Eo 


AEh (jks) 


AEj(jks) 


01s 


±H + 2Jy 


T2H 


+4J X 


lis 


±H + 2{J X + Jy) 


=F2# 





21s 


±H + 2(2J X + Jy) 


T2H 


-47* 


10s f 


±H + 2J X 


=F2# 


+ 4Jy 


20s f 


±H + 4J X 


T2H 


~<i(J X ~ Jy) 


12s X 


±H + 2{J X +2Jy) 


=F2# 


-4J y 


02s \ 


±H + AJy 


T2H 


+ 4(J X - Jy) 


22s * 


±H + 4{J X +Jy) 


=F2# 


-4(Js + Jy) 


00s * 


±H 


T2H 


+ 4(J X + Jy) 
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TABLE II: The mean populations for the spin classes that have nonzero populations in the SOS approximation, with the 
corresponding contributions to the interface velocity under the soft Glauber dynamic. First column: the class labels, jks. 
Second column: the mean spin-class populations for general tilt angle <f>, with cosh 7(0) from Eq. Jjj)- Third and fourth 
columns: the spin-class populations for 4> — [using 7(0) = 0] and <j> = 45° [using Eq. (JSJ for exp[7(45°)]], respectively. Fifth 
column: the contributions to the mean interface velocity in the y direction from spins in classes jk— and jk+, Eq. HUH , using 
the soft Glauber dynamic, Eq. ©. 



Class, jks 


(n(jks)), general <f> 


(n(jks)), = 


(n(jfes)), = 45° 




01s 


1-2X cosii~f(4>)+X 2 


1 


1 


tanh(/3H) 


(1-x 2 ) 2 


(l+xY 2 


2(1+X*) 


l + e 4/JJ x 


lis 


2X[{1+X 2 ) cosh7(»-2X] 


2X 


1 


tanh(/3£T) 




(1+X)2 


2 


2 


21s 


X 2 [1-2X cosh 7(» + X 2 ] 


X 2 


X 2 


tanh(/3ff) 


(1-X 2 ) 2 


(1+X)^ 




l + e -4/3J x 


10s f 


2X 2 f 2cosh 2 7 (0)-l-2X cosh y(<t>) + X 2 


2X 2 (1+2X) 


1+2X 2 +3X 4 


tanh(/3£T) 


1-X 2 1 1-2X cosh 7(0) + X a 

X 2 [1-2X cosh 7(0) + X 2 ] 1 
(1-X 2 ) 2 ] 


(1-X^)(l+Xp 


2(1-X 4 ) 


l + e Wy 


20s f 


X 4 [1-2X cosh 7(^>) + X 2 ] 


X 4 


X 4 


tanh(/3H) 




(1-X^)(1+X)^ 


2(1-X 4 ) 
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FIG. 1: A short segment of an SOS interface y = h(x) between a positively magnetized phase below and a negative phase 
above. The step heights are S(x) — h(x + 1/2) — h(x — 1/2). Interface sites representative of the different SOS spin classes (see 
Table |H and Tabled) are marked with the notation jks explained in the text. Sites in the uniform bulk phases are 00— and 
00+. From Ref. 
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FIG. 2: The average stationary normal interface velocity (v±), shown vs H for <f> = 0. The MC results are shown as data 
points, and the theoretical predictions as solid curves. From below to above, the temperatures are T/T c — 0.2, 0.4, 0.6, 0.8, 
and 1.0. 
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FIG. 3: The average stationary normal interface velocity (v±), shown vs tan0 for (from below to above) H/J — 0.1, 0.5, 1.0, 
1.5, 2.0, 2.5, and 3.0. The MC results are shown as data points, and the theoretical predictions as solid curves, (a) T = 0.2T C . 
(b) T = 0.6T C . 
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FIG. 4: The average stationary normal interface velocity (v±), shown vs T at 4> = for (from below to above) H/J = 0.1, 0.5, 
1.0, 1.5, 2.0, 2.5, and 3.0. The MC results are shown as data points, and the theoretical predictions as solid curves. 
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FIG. 5: Mean stationary class populations {n(jks)), shown vs H at <f> — and T = 0.6T C . Solid lines correspond to the 
theoretical predictions, while filled symbols denote s = +1 and empty symbols (in most cases hidden by the corresponding 
filled symbols) denote s = — 1. (a) The ten SOS-compatible classes, from top to bottom 01s, lis, 10s, 21s, and 20s. (b) The 
six classes with two broken y-bonds, which have zero populations in the SOS approximation, 12s (triangles), 22s (squares), 
and 02s (diamonds), (c) From above to below are shown the aggregate populations of classes with one (01s and 10s), two (lis, 
20s, and 02s), three (21s and 12s) and four (22s) broken bonds. 
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FIG. 6: The two relative skewness parameters p (triangles) and e (squares), denned in Eqs. I|140 and l|15|l . respectively, shown 
vs H for <f> = at T — 0.6T C . (a) The skewness parameters multiplied by 100 for the Ising model with soft Glauber dynamics, 
discussed in this paper, (b) The skewness parameters (no multiplication) for the Ising model with hard Glauber dynamics, 
discussed in Ref. |3|. (c) The skewness parameters (no multiplication) for the SOS model with hard Glauber dynamics, from 



Rcf. 



Note the different sign and the scale difference of two orders of magnitude between the soft dynamic shown in (a) and 



the hard dynamics shown in (b) and (c) 
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FIG. 7: The relative skewness parameter k, shown vs H for <j) — at T — 0.6T C . It measures the asymmetry between the 
populations of bubbles behind and in front of the Ising interface. 
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